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SUMMARY 


This  report  covers  research  and  operations  activities  at 
the  Norwegian  Seismic  Array  (NORSAR)  for  the  period  1 July  - 
30  September  1976.  This  is  the  last  period  where  the  full 
(22  subarray)  array  configuration  is  retained. 

Within  the  reporting  period  there  have  been  170  breaks  in 
the  otherwise  continuous  operation  of  the  Detection  Processor 
(DP).  The  uptime  percentage  is  96.6.  Three  operators  left  in 
the  period,  and  the  vacant  positions  were  covered  by  substi- 
tutes and  overtime  by  the  other  operators.  The  bulletin  pro- 
duction under  the  full-sized  NORSAR  has  now  been  discontinued 
and  some  statistics  for  the  period  1971-1976  are  given. 
Altogether  about  70  000  processings  have  been  completed  by 
the  Event  Processor  (EP)  out  of  which  approximately  36  000 
have  been  accepted  as  true  events.  In  the  DP  a number  of 
modifications  have  been  performed  in  order  to  adapt  the 
system  to  an  environment  with  only  7 subarrays.  At  subarray 
01B  it  has  been  technically  checked  out  that  in  the  new 
array  configuration  it  will  be  possible  to  send  attenuated 
long  period  (LP)  data  in  parallel  to  the  ordinary  data  flow. 
The  performance  of  the  array  instrumentation  has  been  satis- 
factory. Five  reports/papers  have  been  finished  and  four 
topics  are  covered  in  the  summary  of  research  activities.  It 
has  previously  been  reported  on  a project  where  2.778  Hz 
monochromatic  waves  from  a hydroelectric  power  plant  have 
been  used  for  precise  monitoring  of  seismic  velocities.  A 
three-component  short  period  seismometer  set  has  been  in- 
stalled in  site  14C02,  and  using  data  from  these  instruments 
it  is  demonstrated  that  the  waves  have  an  S-type  particle 
motion.  In  the  weighted  beamforming  process  previously 
studied  at  NORSAR,  also  negative  amplitude  weights  were 
allowed.  This  occasionally  led  to  undesirable  side  effects, 

and  in  the  second  study  reported  herein  it  is  described  how 
this  problem  may  be  avoided.  In  the  third  study  simulated 
data  is  used  in  order  to  find  out  what  sort  of  information 


one  may  obtain  about  the  crust  by  using  spectra  of  long 
period  P-waves.  Finally  is  reported  on  a study  of  the 
variance  of  teleseismic  P-wave  amplitudes,  where  it  is 
indicated  that  world-wide  m^  scatter  increases  with  de- 
creasing magnitude,  and  that  magnitude  bias  studies  based  on 
large  events  will  not  necessarily  be  adequate  at  low  magnitudes. 

Three  appendices  are  included  in  this  report.  Appendix  A 
describes  the  format  of  the  NORSAR  High  Rate  tape  while 
Appendix  B describes  the  format  of  the  Low  Rate  tape  for 
data  after  1 October  1976.  In  Appendix  C the  new  array  beam 
deployment  is  described. 


K.-A.  Berteussen 
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OPERATION  OF  ALL  SYSTEMS 


Detection  Processor  Operation 

Within  the  reporting  period  there  have  been  170  breaks  in  the 
otherwise  continuous  operation  of  the  Detection  Processor  (DP) 
system.  The  uptime  percentage  is  96.6%,  as  compared  to  93.8% 
for  the  last  reporting  period  ( January-July  1976) . 

Fig.  II. 1 and  the  accompanying  Table  II.  1.1  both  show  the 
daily  DP  downtime,  in  hours,  for  the  days  between  1 July  and 
30  September  1976.  The  monthly  recording  times  and  up  per- 
centages are  given  in  Table  II. 1.2. 

The  most  significant  break  in  the  recording  occurred  on 
12  September,  when  the  SPS  core  memory  was  overheated  and 
failed.  This  caused  a break  of  more  than  8 hours.  A power 
failure  on  9 August  caused  a stoppage  of  about  3 hours. 

The  170  breaks  occurring  within  the  reporting  period  may 
be  grouped  in  the  following  categories: 


a) 

Error  on  the  multiplexor  channel 

67 

b) 

Software  related  stops 

41 

c) 

SPS  malfunctioning 

17 

d) 

Tests 

12 

e) 

EOC  unit  problems 

8 

f) 

Tape  drive  " 

5 

g) 

TIP  related  stops 

5 

h) 

Customs  Engineering  maintenance 

4 

i) 

Operation 

4 

j) 

Unknown  reason 

4 

k) 

Other  hardware  (TOD,  punch) 

2 

1) 

Power  failure 

1 

4 


The  unrecoverable  error  on  the  multiplexor  channel  (category  a) , 
which  was  described  in  the  report  for  the  last  period,  still 
causes  the  most  frequent  breaks.  So  far  no  solution  has  been 
found  to  this  problem,  which  is  related  to  the  presence  of  the 
TIP  interface  unit  on  this  channel.  However,  tests  performed 
show  that  the  problem  may  be  circumvented  by  writing  the  printer 
output  to  a magnetic  tape. 

As  can  be  seen  from  Table  II. 1.1.  the  software  related  stops 
occurred  frequently  in  the  beginning  of  the  reporting  period, 
to  disappear  towards  the  end  of  the  period.  This  is  related 
to  an  internal  priority  problem  of  the  NCP  task,  which  finally 
was  solved  (see  under  Improvements  and  Modifications,  Chapter 
IV)  . 

The  number  of  SPS  related  stops  is  relatively  normal.  Category 
d)  shows  the  number  of  times  the  system  was  taken  down  to  per- 
form some  kind  of  test  (i.e.,  a new  version,  EOC  hardware  fix, 
etc.).  The  somewhat  large  number  of  tests  was  caused  by  the 
discovery  that  the  Secondary  Online  system  could  no  more  be 
used  simultaneously  with  the  Primary  Online  system,  since 
initialization  of  the  former  invariably  brought  down  the  SPS. 

As  can  be  seen  from  both  tables  (II.  1.1  and  II.  1.2),  the  overall 
situation  was  improved  towards  the  end  of  the  reporting  period, 
when  software  related  stops  were  eliminated.  The  total  down 
time  for  this  period  was  75  hours  39  minutes.  The  mean-time- 
between-f ailures  was  0.5  days  as  compared  with  0.4  days  with 
the  last  reporting  period  (January- July  1976) . 

D.  Rieber-Mohn 
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TABLE  1 1. 1.1 

LIST  CF  8REAKS  IN  CP  PROCESSING  THE  LAST  3 MONTHS 
LAY  START  STOP  COMMENTS 


183 

12 

36 

12 

39 

PROGRAM  HANGLP.RESTART 

18  3 

2C 

21 

20 

37 

MPX/L ATE  ERROR 

183 

22 

43 

22 

50 

MPX/LATl  ERROR 

184 

14 

55 

15 

7 

EOC  RELATED  PROG.  CHECK 

184 

15 

1C 

18 

10 

MPX/LATc,SPS  HARDWARE 

185 

1 

33 

1 

45 

PROGRAM  ERRLR 

185 

10 

42 

10 

48 

EOC  PROBLEMS 

185 

11 

23 

11 

28 

EOC  PROBLEMS 

185 

12 

4C 

13 

13 

MPX/LATL  ERROR 

186 

9 

15 

9 

29 

MPX/LATE  ERROR 

186 

15 

15 

15 

42 

MPX/L ATE  ERROR 

187 

11 

20 

11 

27 

SPS  SYS  HALT 

189 

6 

44 

7 

10 

TAPE  CRIVE  PROBLEM 

189 

8 

4 

8 

17 

MPX/LATE  ERROR 

190 

8 

26 

8 

31 

PROGRAM  ERROR 

190 

11 

54 

12 

•J 

CHNG  OP.  01 SK  PACKS 

190 

14 

20 

14 

24 

PROGRAM  HANGCP 

191 

7 

4 7 

7 

56 

NCP  TASK  HANGUP 

191 

11 

19 

11 

22 

NCP  TASK  HANGUP 

191 

17 

14 

17 

17 

PROGRAM  CHANGE 

191 

18 

59 

19 

4 

NCP  TASK  HANGUP 

192 

13 

26 

13 

34 

NCP  TASK  HANGUP 

192 

13 

50 

14 

11 

TAPE  PROBLEMS 

192 

17 

3u 

17 

53 

MPX/LATE  ERROR 

192 

18 

58 

19 

12 

MPX/LATE  ERROR 

193 

3 

5 

3 

21 

REQUESTED  FROM  CCP 

193 

3 

50 

3 

55 

TIP  PROBLEM 

19  3 

10 

47 

10 

52 

NO  TP  START  TIMES 

193 

20 

32 

20 

43 

TAPE  PRUBLEMS 

194 

5 

52 

6 

6 

SPS  INTER.  NUT  REC. 

195 

0 

12 

0 

29 

SPS 

195 

10 

3 

1C 

fc 

PROGRAM  HANGCP 

195 

12 

19 

12 

29 

PROGRAM  HANGCP 

195 

20 

42 

21 

13 

MPX/LATL  ERROR 

197 

11 

37 

13 

5 

EOC  MAINT  AND  REP. 

197 

13 

58 

14 

15 

CPU  DISABLED 

198 

3 

44 

4 

5 

MPX/LATE  ERROR 

198 

8 

2C 

8 

43 

PRUGRAM  HANGUP 

198 

22 

5 

22 

54 

MPX/LATe  error 

199 

1 

45 

2 

G 

CPU  DISABLEC 

199 

14 

54 

15 

3 

TIP  CRASH 

201 

14 

15 

14 

30 

MPX/LATE  ERROR 

201 

23 

58 

24 

0 

TIP  CRASH 

2C2 

C 

0 

0 

38 

TIP  CRASH 

202 

1 

32 

1 

39 

TIP  CRASH 

202 

11 

48 

11 

54 

PROGRAM  HANGCP 

203 

0 

18 

0 

37 

MPX/LATE  ERROR 

20  3 

10 

8 

10 

53 

CPU  DISABLED 

203 

22 

25 

22 

36 

PROGRAM  HANGCP 

2C4 

9 

2- 

1J 

3 

CNKNOWN  REASON 

204 

10 

50 

11 

9 

CPU  DISABLED 

204 

13 

26 

13 

54 

AM  TEST 

TABLE  II. 1.1 


2C5 

2 

17 

2 

36 

2U5 

15 

26 

15 

33 

205 

17 

51 

17 

58 

206 

7 

24 

7 

57 

2C6 

12 

27 

12 

42 

207 

8 

42 

8 

59 

208 

8 

49 

9 

5 

209 

2 

35 

2 

44 

2C9 

21 

31 

21 

35 

209 

22 

25 

22 

37 

210 

9 

58 

10 

11 

210 

18 

9 

18 

1 / 

211 

1 

16 

1 

31 

211 

3 

4 7 

3 

57 

211 

4 

4 

4 

28 

211 

5 

2C 

5 

35 

211 

1G 

38 

10 

45 

211 

21 

54 

21 

58 

212 

2 

C 

2 

15 

212 

3 

27 

3 

36 

212 

4 

0 

4 

10 

212 

12 

5C 

13 

0 

212 

20 

20 

20 

34 

213 

14 

45 

15 

5 

2 1 A 

1 

6 

1 

1C 

214 

16 

20 

16 

47 

215 

13 

36 

13 

57 

216 

7 

23 

7 

47 

216 

19 

17 

2C 

C 

217 

0 

3 

0 

13 

217 

12 

24 

12 

41 

217 

12 

58 

13 

9 

218 

0 

15 

0 

50 

218 

9 

21 

9 

28 

218 

11 

23 

11 

51 

218 

12 

1 1 

12 

25 

219 

7 

49 

8 

9 

219 

9 

5o 

10 

19 

219 

14 

1 

14 

13 

220 

1 

23 

1 

2e 

220 

11 

45 

11 

54 

221 

8 

12 

8 

25 

221 

9 

48 

10 

2 

221 

16 

35 

16 

46 

221 

18 

5 

18 

20 

222 

9 

15 

9 

46 

222 

11 

59 

14 

37 

223 

13 

52 

14 

1C 

223 

21 

Id 

21 

28 

224 

0 

2 

;> 

W 

9 

224 

10 

2 

10 

20 

224 

14 

24 

14 

40 

224 

16 

0 

16 

40 

MPX/LATE  ERROR 
PROGRAM  HANGLP 
PROGRAM  HANGLP 
MPX/LATE  ERROR 
CPU  DISABLED 
PROGRAM  HANGLP 
AM  TEST 

PROGRAM  HANGLP 
PROGRAM  HANGLP 
MPX/LATE  ERROR 
MPX/LATt  ERROR 
PROGRAM  HANGLP 
MPX/LATE  ERROR 
CPU  DISABLED 
MPX/LATE  ERRLR 
MPX/LATt  ERROR 
PROGRAM  HANGLP 
PROGRAM  HANGLP 
MPX/LATE  ERROR 
TAPE  ORIVE  PROBLEMS 
PROGRAM  HANGLP 
C.  E.  2311 
MPX/LATE  ERROR 
UNKNOWN  REASON 
PROGRAM  HANGLP 
MPX/LATt  ERROR 
MPX/LATE  ERROR 
PRUGRAM  CHECK 
C.  E.  ON  A 2311 
PROGRAM  HANGLP 
CHANGE  OF  DISK  PACK 
EOC  RELATED  STOP 
EOC  RELATED  STOP 
PROGRAM  HANGLP 
MPX/LATE  ERROR 
MPX/LATE  ERROR 
PROGRAM  HANGLP 
CE  ON  1052 
CHANGE  OF  01 SK  PACK 
PROGRAM  HANGLP 
PROGRAM  HANGLP 
TAPE  DRIVE  PROBLEMS 
MPX/LATt  ERROR 
UNKNOWN  REASON 
MPX/LATE  ERROR 
MPX/LATE  ERRUR 
POWER  FAILURE 
SPS  FAILURE 
PROGRAM  HANGUP 
PROGRAM  HANGLP 
PROGRAM  CHANGE 
MPX/LATE  ERRCR 
MPX/LATE  ERROR 
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TABLE  II. 1.1 


227 

2 

3 

2 

18 

MPX/LATt  ERROR 

229 

7 

3 

7 

24 

SPS  F A I LORE 

229 

14 

18 

14 

28 

SPS  FAILURE 

230 

4 

3 3 

5 

16 

MPX/LATt  ERROR 

230 

5 

20 

5 

28 

PROGRAM  FAILURE 

230 

23 

39 

23 

54 

SPS  FAILURE 

232 

6 

19 

6 

51 

MPX/LATE  ERROR 

232 

13 

0 

13 

5 

SPS  FAILURE 

232 

17 

4 

17 

25 

MPX/LATt  ERROR 

232 

22 

41 

22 

58 

MPX/LATE  ERROR 

233 

1 

9 

1 

24 

SPS  FAILURE 

233 

1 

5 . 

2 

11 

SPS  FAILURE 

233 

9 

27 

11 

24 

SPS  ROS  TAPE  CHANGE 

234 

13 

5C 

19 

8 

MPX/LATE  ERROR 

235 

1 

37 

1 

54 

MPX/LATE  ERROR 

236 

10 

38 

1C 

50 

MPX/LATE  ERRCR 

236 

13 

11 

i 3 

28 

CE  ON  1152 

236 

22 

47 

23 

2 

MPX/LATE  ERRCR 

237 

1C 

13 

10 

25 

EOC  MAINTENANCE 

238 

23 

5 1 

24 

G 

MPX/LATE  ERRCR 

239 

0 

0 

2 

MPX/LATE  ERRCR 

239 

23 

58 

24 

0 

MPX/LATE  ERRCR 

240 

0 

0 

0 

16 

MPX/LATE  ERRCR 

24C 

7 

28 

7 

48 

TEST  EOC 

241 

3 

2C 

3 

56 

MPX/LATE  ERRLR 

242 

11 

16 

11 

42 

MPX/LATE  ERROR 

243 

3 

58 

4 

10 

MPX/LATE  ERROR 

243 

14 

15 

14 

23 

SPS  FAILURE 

243 

16 

53 

17 

6 

MPX/LATt  ERRCR 

243 

23 

39 

24 

0 

MPX/LATE  ERRCR 

244 

C 

* 

C 

38 

MPX/LATt  ERROR 

244 

15 

41 

15 

51 

MPX/LATE  ERRCR 

245 

17 

4<„ 

1 / 

58 

UNKNOWN  REASON 

246 

17 

35 

18 

0 

MPX/LATE  ERRLR 

249 

2 

12 

2 

25 

PUNCH  FAILURt 

249 

21 

22 

21 

34 

MPX/LATE  ERRCR 

252 

7 

52 

8 

17 

MPX/LATt  ERRCR 

252 

21 

42 

21 

46 

TOD  FAILURE 

253 

1 

25 

1 

30 

SPS  FAILURE 

255 

15 

50 

lfc 

6 

MPX/LATE  ERRCR 

256 

6 

24 

7 

22 

MPX/LATE  ERRCR 

256 

15 

0 

22 

19 

SPS  CORE  OVERHEATED 

257 

9 

6 

9 

22 

MPX/LATE  ERROR 

257 

17 

4C 

18 

25 

MPX/LATl  ERRLR 

258 

3 

15 

4 

31 

MPX/LATE  ERRCR  t TC . 

259 

13 

16 

13 

34 

SPS  FAILURE 

260 

13 

51 

14 

32 

SPS  FAILURE 

8 


260 

22 

2 

22 

8 

SPS  FAILURE 

261 

14 

6 

15 

40 

TEST  NEW 

VERSION 

263 

13 

30 

14 

34 

MPx/LATE 

ERROR 

264 

12 

16 

13 

42 

TEST  NEW 

VERSION 

26b 

18 

5 

18 

21 

MPX/LATE 

ERROR 

265 

20 

7 

20 

32 

MPx/LATE 

ERROR 

266 

13 

40 

14 

5 

TEST  NEW 

VERSION 

266 

15 

38 

16 

C 

MPX/LATE 

ERROR 

266 

21 

28 

22 

3 

MPX/LATE 

ERRGR 

267 

12 

9 

12 

20 

MPX/LATt 

ERROR 

268 

9 

51 

10 

34 

TEST  NEW 

VERSION 

268 

11 

55 

13 

11 

TEST  NEW 

VERSION 

268 

23 

54 

24 

0 

MPX/LATE 

ERROR 

269 

C 

C 

C 

27 

MPX/LATE 

ERROR 

269 

7 

45 

8 

31 

MPX/LATE 

ERROR 

270 

7 

23 

8 

23 

MPX/LATE 

ERROR 

270 

13 

27 

13 

31 

TAKEN  DOWN  8Y  OPERATOR 

271 

8 

14 

9 

3 

TEST  NEW 

VERSION 

271 

12 

52 

13 

34 

TEST  NEW 

VERSION 

271 

18 

14 

18 

46 

MPX/LATE 

ERROR 

272 

9 

11 

10 

47 

TEST  NEW 

VERSION 

273 

8 

1 

iu 

57 

TEST  NEW 

VERSION 

274 

11 

34 

12 

9 

START  NEW 

i VERSION 

I 


DP  and  EP  computer  usage  July-September  1976 
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11 . 2 Event  Processor  Operation  (EP) 

The  Event  Processor  system  has  performed  satisfactorily 
throughout  the  reporting  period.  Its  up  time  percentage 
is  24.3,  as  compared  to  29.6  for  the  last  reporting  period 
( January-June  1976) . 

H.  Bungum 

11 . 3 NORSAR  Data  Processing  Center  (NDPC)  Operation 
Data_Center 

The  data  center  operation  has  been  affected  by  the  planned 
change  to  unmanned  operation.  Three  operators  left  in  the 
period.  The  vacant  positions  were  covered  by  substitutes 
and  overtime  by  the  other  operators. 

The  Detection  Processor  up  time  was  96.6%  for  the  period. 
Here  bugs  in  and  testing  of  the  new  DP  version  stand  for 
about  1%  of  the  down  time. 

The  SPS  had  two  major  stops,  both  times  it  was  put  back  into 
operation  by  project  personnel.  Also  on  the  EOC  project  per- 
sonnel have  been  engaged  in  f ault-f inding  maintenance. 

The  power  situation  was  fairly  good,*  thunderstorms  usually 
cause  a number  of  power  breaks  in  this  period. 

2§ta_Communication 

The  Terminal  Interface  Message  Processor  (TIP) 

After  a power  outage  9 August  the  TIP  could  not  be  restarted 
and  Mr.  Kelly  (Bolt,  Beranek  and  Newman)  had  to  be  called 
in  from  London.  After  this  incident  the  TIP  has  performed 
as  normal.  The  power  supply  was  replaced  23  August.  Further 
corrections  to  the  TIP  have  not  been  made  since  the  last 
report. 
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Na  t ional_Communicafc  ion_C±r;Guits 

Outages  have  occurred  sporadically  and  groups  of  subarrays  have 
been  affected.  Individual  subarrays  have  also  been  subject  to 
outages  or  deteriorated  performance  due  to: 

Cable  faults  (06B) 

Line  trouble/EPU  power  interference  (01C) 

Broken  cable  (02C) 

NTA-work  (PCM-installation  in  the  area)  ( 02B, 03C, 05C) 
Intermittent  operating  equal izer/ampl . (06C) 

Broken  cable/lack  of  power/faulty  equalizer/amplifier 
(IOC) 

lD^SED3tional_Communication_Circuits__(Londoni_SDAC)_ 

The  London  CCT 

Apart  from  minor  irregularities,  there  was  a reliable  per- 
formance this  period. 

The  SDAC  Circuit 

The  performance  of  this  circuit  has  been  changing,  and  has 
been  subject  to  tests  quite  frequently.  These  tests  have  been 
initiated  by  NCC  (Network  Control  Center)  but  also  by  ITT 
July  was  a quite  good  period  with  few  irregularities.  For 
August,  however,  the  picture  changed:  4,  5,  9,  10  and  20 
August  were  affected  by: 

Modem  line  test  (NCC  request) 

Loss  of  carrier 

Temporary  line  rerouting  (NTA/Lillestr0m) 

Loss  of  carrier  and  frequently  modem  tests  in 
different  modes  (NCC  request) 

Modem  check  (4  hours) , respectively. 

The  last  month,  September,  has  the  following  incidents: 

1 September  different  line  tests  were  requested  by  ITT  after 
SDAC  had  claimed  errors  in  data  exchange  with  NORSAR.  NTA 
and  NDPC  personnel  were  involved.  According  to  indicators 
on  the  SDAC  modem  the  line  parameters  were  not  outside  speci- 
fications. 8 September  the  SDAC  modem  was  put  into  "Audio 


Bus  Back"  after  NTA  request,  i.e.,  line  loop  in  front  of  the 
modem  at  NDPC.  14  September  message  from  NCC  that  line  checks 
would  be  carried  out  in  connection  with  introduction  of  a new 
"16  channel  system",  which,  if  successful,  will  improve  the 
system's  capability  to  handle  more  "packets"  (messages). 


TABLE  I I. 3.1 


Communications,  degraded  performance  ( >20) /outages  (>200).  Figures  in 
per  cent  of  total  times.  Month  = 4 or  5 periods  as  indicated. 


August (4) 

I >200 


(03B) 


1.0 

(10C) 

II .4  ARPANET 

There  has  been  no  change  in  the  attachment  configuration  to, 
or  use  of,  the  NORSAR  TIP  since  the  last  reporting  period. 


D.  Rieber-Mohn 


ARRAY  PERFORMANCE 


Table  III.l  shows  some  basic  statistics  for  the  EP  operation 
during  the  present  reporting  period,  in  terms  of  analyst  de- 
cisions for  all  the  DP  detections  reported  by  EP.  Table  III. 2 
shows  the  number  of  reported  events  on  a monthly  basis,  and 
in  Fig.  III.l  is  shown  the  distribution  on  a daily  basis. 

The  numbers  are  quite  normal,  with  an  average  of  20.7  events 
per  day. 

Complete  Statistics  1971-76 

The  bulletin  production  under  the  full-sized  NORSAR  has 
now  been  discontinued,  and  it  could  therefore  be  interesting 
to  look  at  the  complete  statistics  up  to  30  September  1976. 

I 

The  EP  processing  statistics  are  given  in  Table  III. 3,  where 
data  are  available  from  1 July  1971.  It  is  seen  there  that 
about  72  000  processings  have  been  completed  by  EP,  out  of 
which  50%,  or  36  000,  have  been  accepted  as  true  events.  Note 
that  the  processing  time  was  systematically  decreased  to 
a final  value  of  about  15  minutes  per  event,  with  a correspond- 
ing  decrease  in  the  EP  processing  load  on  the  B-computer . This 
increase  in  efficiency  was  the  result  of  a systematic  and  con- 
centrated effort  on  the  part  of  the  NTNF/NORSAR  personnel  in 
order  to  improve  the  processing  packages  and  procedures  originally 
developed  and  delivered  by  IBM/Federal  Systems  Division. 

In  Table  III.  4 we  have  given  the  number  of  reported  events 
(teleseismic  and  core  phases)  for  each  year  from  1 May  1971 
to  30  September  1976.  The  reason  why  the  total  number  (33809) 
is  somewhat  smaller  than  that  given  in  Table  III. 3 for  the 
number  of  accepted  processings  (35758)  is  that  the  latter 
sometimes  includes  several  processed  phases  for  the  same  event. 

It  is  seen  that  the  overall  average  per  day  is  17.1,  for  the 

data  after  1973  that  number  is  18.5. 


H.  Bungum 


I 
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Analyst 

Classification 

No.  of 
Processings 

Percentage 

Accepted  as  events 

2070 

47.8 

Rejected  as  being: 

- Poor  SNR  or  noise 

447 

10.3 

- Local  events 

713 

16.5 

- Double  jrocessings 

559 

12.9 

- Communication  errors 

542 

12.5 

Sum  Processed 

4331 

100.0 

l 

TABLE  III. 2 

Number  of  teleseismic  and  core  phase  events  reported  during  the  time 

period  July-Sept  1976 
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TABLE  1 1 1 . 3 

Statistics  of  the  total  EP  history  from  1 July  1071  to  30  Sept  1076 


Time 

Period 

EP% 

— 

EP 

Hours 

No . of 
Proc. 

Time  per 
Proc.  (min) 

No.  of 
Acc.  Proc. 

Acceptance 

Percentage 

Time  per 
acc . event 
(min) 

1971-2 

37.9 

1673 

1706 

58.8 

971 

56.9 

103.4 

1972-1 

39.7 

1736 

4574 

22.8 

2913 

63.7 

35.8 

-2 

45.8 

2023 

8070 

15.0 

3729 

46.2 

32.6 

1973-1 

40.3 

1749 

7802 

13.5 

3899 

50.0 

26.9 

-2 

41.1 

1815 

8176 

L 3 . 3 

3478 

42.5 

31.3 

1974-1 

29.5 

1283 

7023 

11.0 

3298 

47.0 

23.3 

-2 

32.5 

1434 

7845 

11.0 

3798 

48.4 

22.7 

1975-1 

27.2 

1175 

7380 

9.6 

3929 

53.2 

17.9 

-2 

25.5 

1128 

6436 

10.5 

3214 

49.9 

21.1 

I 1976-1 

29.6 

1292 

84  59 

9.2 

4459 

52.7 

17.4 

-2 

24.3 

536 

4331 

7.4 

2070 

47.8 

15.5 

j 

Sum 

15844 

35758 

. Daily 

34.4 

8.3 

37.4 

13.2 

18.6 

26.6 

Average 



TABLE  I I I. 4 


Number  of  events  reported  by  NORSAR  on  a yearly  basis,  from 
1 May  1971  to  30  Sept  1976 
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IV. 
IV.  l 


IMPROVEMENTS  AND  MODIFICATIONS 
Detection  Processor 

Within  this  period,  further  modifications  have  been  imple- 
mented in  the  Detection  Processor  system.  These  modifications 
fall  into  two  categories: 

a.  general  error  corrections  and  improvements, 

b.  modifications  to  adapt  the  system  to  an  environment  with 
7 subarrays  only  (versus  22  earlier)  . 

Under  a)  the  following  changes  were  done  to  the  system. 

The  intra-  and  inter-task  flow  of  ARPANET  messages  was 
considerably  simplified.  This  was  achieved  by  using  one 
and  the  same  queue  block  for  a message  for  the  duration 
of  its  flow  through  the  system.  In  this  way,  the  over- 
head of  releasing  and  leasing  queue  blocks  was  reduced, 
and  the  algorithms  became  simpler. 

The  handling  of  timed-out  messages,  waiting  for  acknowledge- 
ment (RFNM) , was  given  to  a separate  subtask,  which  will 
be  activated  periodically,  with  intervals  equal  to  the 
timeout  period  for  messages  sent  out.  Any  stop  in  the 
flow  of  outgoing  messages  because  of  non-removal  of  timed- 
out  messages  should  thus  not  occur. 

The  problem  of  giving  (internal)  priorities  to  the  various 
subtasks  within  the  NCP  task  was  finall'  solved,  by  adopt- 
ing the  principle  of  giving  highest  priority  to  subtasks 
with  the  lowest  frequency  of  execution.  In  this  way,  high 
priority  was  given  to  the  subtask  handling  timed-out 
messages  (IMPTIM) , the  task-initial ization-and-operator- 
message  subtask  (SNNCPR) , and  tne  Imp-to-Host  command 
interpreter  subtask  (IMPHOS),  while  the  subtasks  handling 
the  in-  and  outgoing  ARPANET  messages  (IMPREAD,  IMPIN, 

IMPOUT  and  IMPWRITE)  got  low  priority.  This  setup  seems 
to  prevent  any  deadlock  in  the  flow  of  ARPANET  messages 
within  the  task. 
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Numerous  modifications  have  been  made  in  other  tasks, 
as  queuing  Incoherent  Detections  for  ARPANET  transmission, 
verification  of  tape  start/stop  times,  removal  of  all  code 
related  to  the  EOC  Beam  Display  (not  used  any  more)  and 
removal  of  all  reference  to  the  EOC  debug  tape. 

Under  b)  the  following  changes  and  additions  were  made  to  the 
system. 

The  PNRSPS  processor,  which  reads  the  SPS  block,  was 
modified  to  read  High  Rate  data  from  14  subarrays  only 
(this  because  06C  - the  last  subarray  in  the  new  con- 
figuration - is  the  14th  subarray  in  the  old  setup) . 

The  PNRNOR  processor  was  modified  to  insert  shortened 
'manual'  status  fields  from  DP  Global  Common.  Also,  tests 
on  instruments  for  DC  offset,  etc.,  will  be  performed 
on  valid  subarrays  only. 

The  PNRSAD  processor,  which  sets  up  data  messages  for 
ARPANET  transmission,  was  changed  to  extract  status  and 
data  for  the  7 valid  subarrays  only,  for  insertion  into 
the  transmission  message  block. 

The  PNLRTP  processor,  which  produces  the  Low  Rate  tape, 
was  modified  to  call  the  subroutine  LRXTRT,  which  extracts 
LP  data  from  the  High  Rate  block  for  the  7 valid  subarrays, 
and  builds  a blocked  Low  Rate  NR-record. 

The  PNHRTP  processor,  which  produces  the  High  Rate  tape 
was  modified  to  call  the  COMPRESS  subroutine,  which 
extracts  SP  data  from  the  High  Rate  block  for  the  7 valid 
subarrays  and  builds  a blocked  High  Rate  'NR '-record. 

This  record  will  be  written  to  tape  when  filled. 

A new  subroutine  LRXTRT  (coded  in  PL360)  was  added  to 
the  tape  task.  This  subroutine  extracts  LP  status  and 
data  from  integer  second  High  Rate  blocks,  builds  a Low 
Rate  logical  record,  inserts  this  record  into  a physical 
record  buffer  and  writes  the  buffer  to  tape  when  it  is 
full  and  a write  operation  is  possible. 


- 21  - 


A new  subroutine  COMPRESS  (coded  in  PL360)  was  added  to 
the  tape  task.  This  subroutine  extracts  SP  status  and 
data  from  the  current  HR-block  received  from  the  Data 
Acquisition  task,  builds  a compressed  logical  record  and 
inserts  this  into  the  buffer  for  the  High  Rate  'NR' 
record . 

The  PNRERR  processor  was  modified  to  print  the  periodic 
Error  Reports  for  valid  subarrays  only,  and  only  if 
any  of  the  error  counts  for  the  passed  interval  were 
non-zero . 

The  size  of  the  largest  queue  block,  defined  in  RMONCOMS, 
was  reduced  from  1250  to  775  FWs . This  could  be  done  because 
the  latter  size  was  adequate  for  the  transportation  of 
HR-data  from  14  SAs  through  the  system  (instead  of  22). 

A separate,  fixed  buffer  capable  of  storing  up  to  4 com- 
pressed logical  HR  records  before  the  write  operation, 
was  added  to  the  system  (HRBUFFS) . 

D.  Rieber-Mohn 


IV . 2 Event  Processor 

No  modifications  were  done  to  the  EP  system  during  this  period. 


D.  Rieber-Mohn 


Array  Instrumentation 

The  project  described  in  the  last  semiannual  report  (1  January 
30  June  1976)  of  running  a three-axis  SP  seismometer  at 
subarray  14C  has  continued  with  good  results,  except  for  noise 
on  the  Z-component  data,  LF  oscillations  in  the  period  24  July 
to  5 August  due  to  water  in  the  seismometer  plug  and  super- 
imposed 50  Hz  noise,  corrected  24  August.  As  previously, 
channel  03  and  the  LP  channels  at  14C  have  been  attenuated 
30  dB,  and  channel  14C05  has  been  used  for  phase  shift 
reference  (Section  VII. 1,  this  report).  The  lowpass  active 
filter  at  the  LTAs  were  bypassed  on  the  channels  14C01,  02, 

04  and  06  in  this  period. 

The  array  reduction  program  started  in  September  with  cutting 
and  plugging  of  the  CTV  air  pipes  on  five  sites. 

In  the  new  array  configuration  attenuated  LP  channels  are 
planned.  If  a sampling  rate  of  every  two  seconds  are  suffi- 
cient, it  is  possible  to  sample  the  output  from  the  three 
LP  seismometers  at  a site  in  parallel , by  feeding  the  (attenu- 
ated) signal  to  random  data  addresses  (RDA)  13,  14  and  15  for 
vertical,  NS  and  EW  respectively.  It  was  tecnically  checked 
out  at  01B  from  15  September,  by  feeding  NS  seismometer  data 
output  to  RDA13 . 


Alf  Kr . Nilsen 
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V . 


MAINTENANCE  ACTIVITY 


This  section  includes  a review  of  the  maintenance  accomplished 
at  the  subarrays  by  the  field  technicians  as  a result  of  the 
remote  array  monitoring  and  visual  inspections.  There  are  no 
changes  in  the  monitoring  schedule  this  period,  but  due  to  fre- 
quent EOC  and  DP  faults  the  schedule  has  been  only  loosely 
followed . 


Maintenance  Visits 

Fig.  V.l  shows  the  number  of  visits  to  the  subarrays  in  the 
period.  Excluding  visits  caused  by  troubles  in  the  communi- 
cation system,  the  subarrays  have  on  the  average  been  visited 
2.8  times  (excluding  OlA  and  06C:  2.1  times).  The  large  number 
of  visits  to  OlA  are  due  to  painting  of  the  LPV  and  cable 
breakages,  at  06C  due  to  troubles  with  the  SLEM  and  cable 
breakage . 


No . of 
Visits 


OlA  02B  04B  06B  01C  03C  05C  07C  09C  11C  13C 

01B  03B  05B  07B  02C  04C  06C  08C  10C  12C  14C 

. Fi9-  v-i  Number  of  maintenance  visits  to  the  NORSAR  subarrays 

1 July  - 30  September  1976. 


Preventive  Maintenance  Projects 

The  preventive  maintenance  work  in  the  array  is  described  in 
Table  V.l. 


TABLE  V.l 


1 

j Unit 
— 

Action 

No . of 
Accomp. 

Actions 

Remaining 

Comments 

j LTA 

Adjustment  of  SP  DC 
offset 

11 

Adjustment  of 
channel  gain  SP 

7 

Adjustment  of 
channel  gain  LP 

1 

. 

Power 

Battery  Maintenance 

3 

LPV 

Painting  of  LPV* 

0 

, 1 

01A , 02C 

I 

!*  Reported  as  corrective  maintenance  in  the  montly  reports  due  to 

the  great  need  for  this  work,  but  have  preventive  effects  as  well. 
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Disclosed  Malfunctions  on  Instrumentation  and  Electronics 

Table  V.2  gives  the  number  of  accomplished  adjustments  and 
replacements  of  field  equipment  in  the  total  array  with  the 
exception  of  those  mentioned  in  Table  V.l. 


table  v.2 

Total  number  of  required  adjustments  and  replacements  in  the  NORSAR 
data  channels  and  SLEM  electronics 
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Malfunction  of  Rectifiers,  Power  Loss,  Cable  Breakages 

Malfunction  of  the  rectifiers  and  power  loss  requiring  action 
of  the  field  technicians  or  local  power  company  are  reported 
in  Table  V.3. 


TABLE  V. 3 

Faults  disclosed  in  subarray  rectifiers  and  power  loss 


Sub- 

Array 

Fault 

Period  of 
Inoperation 

Comments 

01C 

No  power 

_ II  _ 

at  site 

8-9  July 
12-12  July 
15-31  July 

Rectifier  did  not  re- 
start after  AC  power 
break . EPU  converter 
card  faulty. 

06C 

Noise  on 

data 

10-22  Sept. 

EPU  and  DU  replaced 

10C 

No  power 

at  site 

12-13  Aug 

Blown  main  fuses  and 
one  fuse  in  battery 
charger 

12C 

No  power 

to  MP  lamps 

4v01t.ronic  power  sup- 
ply replaced. 
— 

Nine  cable  breakages  have  been  repaired,  or  attempts  to  repair 
have  been  made,  requiring  17  days'  work  of  the  field  technicians. 
Due  to  the  short  time  left  before  reduction  of  the  array, 
the  repairs  were  not  completed  at  04B02  and  07C03.  At  01A04 
the  cable  repair  was  dropped,  as  we  intend  to  use  this  channel 
for  attenuated  SP  data. 

Conclusion 

The  performance  of  the  array  instrumentation  has  been  satis- 
factory in  the  period  and  compared  with  previous  periods  no 
anomalies  in  the  operation  were  observed. 


- 27  - 


Two  power  supplies  of  NORSAR  analog  SP  station  (NMC/04B05) 
were  damaged  by  lightning  10  July.  Difficulties  in  acquiring 
spare  parts  have  caused  inoperation  from  that  date  and  out 
the  period. 

A.  Kr.  Nilsen 


ABBREVIATIONS 


AC 

- 

Alternate  current 

ADC 

- 

Analog-to-Digital  Converter 

BB 

- 

Broad  band 

BE  card 

- 

Lightning  protection  card 

CMR 

- 

Common  mode  rejection 

CTV 

- 

Central  Terminal  Vault 

DC 

- 

Direct  current 

DP 

- 

Detection  processor 

DCO 

- 

DC  offset 

DU 

- 

Digital  unit 

EOC 

- 

Experimental  Operations  Console 

EPU 

- 

External  power  unit 

LF 

- 

Low  frequency 

LP 

- 

Long  period 

LPV 

- 

Long  period  vault 

MP 

- 

Mass  position 

NMC 

- 

NORSAR  Maintenance  Center 

RA-5 

- 

SP  seismometer  amplifier 

RCD 

- 

Remote  centering  device 

RSA 

- 

Range  switching  amplifier 

SLEM 

- 

Seismic  short  and  long  period  electronics  module 

SP 

- 

Short  period 

WHV 

- 

Well  head  vault 
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L.B.  Tronrud 


VI. 2 Program  Documentation 


No  program  documentation  has  been  written  during  this  period 


D.  Rieber-Mohn 
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VII.  SUMMARY  OF  SPECIAL  TECHNICAL  REPORTS/PAPERS  PREPARED 

VI I. 1 Seismic  Velocity  Variations  and  Solid  Earth  Tides 

We  have  previously  (Final  Technical  Report,  1 July  1974- 

30  June  1975 ",  Semiannual  Technical  Report  1 July  - 31  December 

1975)  reported  on  this  project,  in  which  2.778  Hz  monochromatic 

waves  from  a hydroelectric  power  plant  have  been  used  for 

precise  monitoring  of  seismic  velocities.  It  has  been  found 

that  precisions  of  the  order  of  10  ^ are  obtained  when  2 

-4 

hours  of  data  are  stacked,  and  10  with  7 days  of  data.  The 
average  group  velocity  is  around  3.5  km/s. 

In  order  to  determine  the  wave  type  that  we  are  observing, 
a three- component  set  of  short  period  seismometers  was  in- 
stalled in  site  14C02,  4.7  km  away  from  the  source.  The  particle 
motion  was  then  derived  by  application  of  Fourier  analysis 
and  stacking,  where  relative  amplitudes  and  phase  delays  were 
used  in  the  actual  contructions  of  the  plot,  as  shown  in 
Fig.  VII. 1.1.  It  is  seen  there  that  the  radial  to  transverse 
amplitude  ratio  is  so  small  that  only  S-waves  can  satisfactorily 
explain  the  observations.  This  solution  is  also  consistent  with 
the  velocity  observations  around  3.5  km/s. 

A systematic  search  for  periodicities  in  the  data  has  revealed 
the  presence  of  velocity  variations  with  a cycle  time  of 
12  hours.  This  is  shown  in  Fig.  VII. 1.2,  where  the  two  full 
cycles  in  the  24  hours  plot  demonstrate  tii^..  we  are  not 
observing  a hidden  diurnal  component  (such  as  a non-symmetric 
temperature  effect) . We  have  searched  the  data  for  a number 
of  other  periods  too,  and  find  a response  for  periods  in  the 
range  of  11-13  hours.  This  of  course  should  be  expected  when 
only  7 days  of  data  have  been  analyzed,  but  it  is  noteworthy 
here  that  the  measured  phase  delays  for  the  different  periods 
in  the  11-13  hours  range  are  found  from  simulation  analysis 
not  to  be  consistent  with  leakage  from  a single  spectral  com- 
ponent . 
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The  closest  place  to  look  for  a source  of  these  velocity 
variations  is  in  the  earth  tidal  field,  which  has  been 
calculated  theoretically  through  the  harmonic  expansion  of 
the  tide-generating  potential  of  second  order.  In  this  way, 
the  horizontal  components  of  the  strain  tensor  have  been 
calculated  for  the  time  periods  in  question,  but  no  correla- 
tion has  been  found  with  the  observed  velocity  variations, 
with  the  important  exception  that  the  periods  coincide.  This 
lack  of  correlation  so  far  as  initial  phase  is  concerned  does 
not,  however,  exclude  the  possibility  of  a tidal  source.  This 
is  because  there  are  factors  such  as  ocean  loading  and  at- 
mospheric pressure  variations  which  have  not  been  taken  into 
consideration  in  the  theoretical  calculations.  More  important, 
however  is  the  doubt  that  we  have  as  to  the  assumption  of 
a direct  relationship  between  velocity  and  strain,  this  is 
partly  because  the  (unknown)  distribution  in  direction  and 
size  of  the  cracks  in  the  material  also  enters  this  relationship. 

The  velocity  variations  displayed  in  Fig.  VII. 1.2  are  of  the 
order  of  10  . The  predictions  here,  based  on  extrapolation 

of  laboratory  values,  give  values  of  the  order  of  10  ^ . If 
the  velocity  variations  that  we  observe  are  actually  caused 
by  earth  tides,  this  therefore  means  that  there  is  a very  high 
velocity  sensitivity  to  stress  variations  in  the  uppermost 
layer  of  the  earth. 


H . Bungum 
T.  Risbo 
E.  Hjortenberg 
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Fig.  VII. 1.2  Propagational  phase  angle  differences  when  measured 
with  a sliding  window  procedure  in  such  a way  that 
the  presence  of  the  searched-for  peiiodicity  should 
show  up  in  the  form  of  a one  cycle  sine  wave.  This  is 
the  case  for  T=12  hours,  while  it  is  the  12  hours 
periodicity  which  shows  up  as  two  cycles  on  the  plot 
for  T=24  hours. 
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VII. 2 Weighted  Beamf orming  using  Non-negative  Weights 

Considerable  attention  has  been  given  to  the  weighted  beam- 
forming process  of  Christoff ersson  and  ilusebye  (1974),  where 
negative  weights  were  allowed.  Fyen  et  al  (1975)  showed  that 
the  presence  of  negative  weights  permitted  additional  statis- 
tical tests  for  use  by  the  analyst  in  the  signal-noise  classi- 
fication. The  weights  used  in  this  beamforming  process  arose 
from  the  model 

Yj  = YjS  + nj  j = 1, . . . ,M  (1) 

where  y^  is  the  recorded  data  at  sensor  j,  S is  the  unknown 
signal  and  Yj  an  unknown  scaling  factor  which  accounts  for  the 
observed  amplitude  variation  at  NORSAR.  n^  is  the  residual 
noise  of  the  model.  The  least  squares  estimators  of  this  model 
are 


s = 1$  j 


(2) 


where  y^  are  the  elements  of  the  eigenvector  corresponding 
to  the  largest  eigenvalue  of  the  matrix  y'y,  i.e.,  the 
varianc 2-covariance  matrix  of  the  recordings  (first  principal 
component) . 


When  ^ takes  negative  values,  it  means  physically  that  the 
sensor  output  has  a phase  shift  of  it  relative  to  the  beam. 

This  in  turn  could  occasionally  lead  to  undesirable  side  effects 
as  the  phase  shift  option  may  result  in  projection  of  pure  noise 
wavelets  into  the  signal  space.  A project  was  therefore  carried 
out  where  the  estimated  weights  were  permitted  to  take  non- 
negative values  only.  The  method  of  least  squares  leads  to 
minimizing 
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M M 

F = F (Y.*  - YiS)'(yi  - Y,S)  - A(  y y-  - D (3) 

j=l  3 3 33  j=l  3 

wjjere  X is  the  Lagrangian  multiplier  reflecting  the  restriction 
minimum  of  F is  found  by  an  iteration  procedure 
0 is  permitted.  This  gives  the  same  solution 
as  the  principal  component  method  when  all  y^  > 0 are  positive. 
However,  when  some  of  the  y-  took  negative  values,  the  iteration 
procedure  will  give  a solution  which  gives  zero  weights  to 
sensors  out  of  phase  and  slightly  different  values  for  others. 

The  procedure  was  applied  to  a set  of  mostly  weak  earthquakes 
with  varying  signal-to-noise  ratios  and  extreme  noise  cases 
which  trigger  the  on-line  detector  (false  alarms) . The  results 
showed  that  non-negative  weights  do  not  give  an  opportunity 
to  create  additional  statistical  tests  for  signal-noise  classi- 
fication. However,  when  calculating  signal-to-noise  (SNR)  ratios, 
the  analysis  showed  that  non-negative  weights  preserve  the 
same  SNR  as  the  principal  component  procedure,  and  for  extreme 
noise  cases  the  SNR  is  reduced  by  at  least  0.15  units  on  an 
average.  In  this  region  of  SNR  we  know  that  the  number  of  detec- 
tions increases  exponentially  with  decreasing  SNR,  so  a reduction 
of  SNR  for  extreme  noise  of  0.15  units  is  considered  as  a signi- 
ficant improvement  relative  to  principal  component  beams,  in 
terms  of  an  SNR  detector. 
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J . Fyen 

A.  Christoff ersson,  Uppsala 
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I . 3 Long  Period  P-Wave  Spectra  as  a Tool  for  Studies  of  Local 
Structure 

The  vertical  (Z)  and  horizontal-radial  (R)  signal  spectra  for 
a body  wave  may  be  written 


Z = S • Hz  • I + Nz 


(1) 


R = S • Hr  • I + Nr 


(2) 


where  S is  the  source  and  the  source-side  effects,  H„ (H„) 
is  vertical  (horizontal)  transfer  function  of  the  crust  and 
upper  mantle  on  the  receiver  side,  I is  instrument  response 
and  Nz  (NR)  is  noise  on  the  vertical  (horizontal)  component. 
All  the  above  parameters  are  a function  of  the  angular  fre- 
quency. Assuming  N = N = 0 one  may  form  the  ratio 

G = Z/R  = Hz/Hr  (3) 


Phinney  (1964)  originally  proposed  the  idea  of  dividing  the 
spectra  of  the  vertical  component  of  a recorded  P-wave  with  that 
of  the  horizontal  radial  component  in  this  way,  thereby  obtaining 
a ratio  which  does  not  depend  on  the  spectrum  of  the  original 
incident  pulse.  These  ratios  are  then  compared  with  theoretical 
layered  models  using  the  Haskell  ( 1953 ) -Thomson  (1950)  matrix 
theory.  This  'spectral  ratio'  method  has  sin  e been  used  by  a 
large  number  of  authors. 


The  aim  of  this  work  (Berteussen,  1976a)  has  been  to  find  out 
what  type  of  information  one  may  expect  to  get  about  the  struc- 
ture under  a station  (having  a three-component  long  period 
instrument  set)  by  applying  this  spectral  ratio  method.  r"o  do 
this  we  have  used  only  simulated  data,  while  we  in  a later 
work  (Berteussen,  1976b)  also  will  try  real  data  from  the 

three-component  long  period  instruments  at  the  Norwegian 
Seismic  Array  (NORSAR) . 
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Simulating  a number  of  different  models,  we  have  found  that  by 
far  the  most  pronounced  effect  on  the  spectral  ratio  is 
caused  by  the  Moho  boundary.  For  this  boundary  the  effect  of, 
for  example,  changing  the  depth  is  very  obvious  and  can  easily 
be  seen  on  the  simulated  ratios.  Other  interfaces  may  compli- 
cate the  picture,  like  for  example  enlarging  one  of  the  peaks 
caused  by  the  Moho  boundary,  but  it  is  extremely  difficult 
without  other  evidence  to  assign  these  additional  effects 
to  a particular  model  (Fig.  VII. 3.1  and  VII. 3. 2).  For  example, 
interfaces  at  10  km  depth  and  30  km  depth  can  give  rise  to 
almost  the  same  set  of  peaks.  Compare  trace  4 from  top  on 
Fig.  VII. 3.1)  with  traces  4 and  6 on  Fig.  VII. 3. 2.  Structures 
below  Moho  may  give  a large  number  of  additional  smaller 
peaks  (top  traces  Fig.  VII. 3. 4),  but  neither  of  these  structures 
seem  able  to  change  drastically  the  spectral  ratio  generated 
by  the  Moho  boundary.  Also  it  may  be  difficult  even  for  a one- 
layer  case  to  discriminate  between  a deep  crust  with  high 
velocity  and  a shallow  one  with  lower  velocity  (Fig.  VII. 3. 3). 

By  taking  into  account  also  the  size  of  the  spectral  peaks, 
this  should,  however,  be  possible  although  the  peak  sizes  are 
easily  affected  by  noise. 

When  using  Fourier  theory  to  estimate  the  vertical  and  hori- 
zontal spectra,  we  will  have  to  use  a time  window  of  a certain 

length.  The  estimated  ratio  is  then  (a  star  means  convolution) 

« _ Z*W  _ (S'VI+NZ>*W 

U R*W  (S-Hd*I+Nd) *W  (4' 

K K 

Since  convolution  is  an  integral  operation,  we  can  in  general 
not  cancel  the  common  terms  in  the  numerator  and  denominator 
in  Eq.  (4),  even  if  we  assume  Nz=NR=0.  In  the  practical  esti- 
mation of  the  spectral  ratios,  one  has  to  ensure  that  a long 
enough  window  is  used  in  the  calculation  of  the  spectra  (see 
Fig.  VII.3.4).  For  the  type  of  instruments  used  at  NORSAR  this 
means  50  seconds  or  longer,  but  with  other  instruments  with  a 
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flatter  response  curve  a somewhat  shorter  window  may  be  used, 
although  as  demonstrated  a 20  second  window  is  surely  too  short. 
If  one  would  like  to  see  the  full  effect  of  the  upper  mantle 
on  these  ratios,  a window  of  up  to  180  sec  is  needed.  There 
are  several  other  factors  to  be  aware  of.  One  needs  signals 
with  very  good  signal-to-noise  ratio,  and  the  signal  should 
be  as  pulse-like  as  possible.  A very  deep  source  is  to  be 
preferred  (one  thereby  may  avoid  any  source-side  depth 
phases),  but  is  not  a completely  necessary  requirement.  It  is 
found  that  depth  phases  are  only  a problem  when  they  are  so 
strong  and  arriving  so  late  that  they  may  be  seen  on  the  record. 
By  selecting  only  simple  pulse-like  signals  one  therefore 
can  avoid  this  problem. 

In  short,  from  the  simulations  performed  we  conclude  that  one 
can  use  the  spectral  ratio  method  to  find  the  depth  to  Moho 
if  proper  care  is  taken  in  the  data  processing.  We  do  not 
believe  that  it  is  possible  to  find  a more  detailed  model, 
at  least  not  without  any  other  evidence. 

K.-A.  Berteussen 
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Fig.  VII. 3.1  Illustration  of  the  effect  on  the  theoretical  spectral  ratio 
of  low  velocity  surface  layers,  in  the  first  trace  is  plotted 
the  ratio  for  a one-layered  reference  model  with  crustal  P- 
velocity  equal  to  6.G  km/sec  and  Moho  depth  35  km.  In  the 
following  cases  are  then  plotted  the  ratios  for  different 
cases  of  surface  layers.  The  thickness  of  the  layer  and  its 
P-velocity  (in  brackets)  is  written  on  the  figures. 


Illustration  of  the  effect  on  interfaces  in  the  deeper  part 
of  the  crust.  As  before  the  upper  trace  is  the  spectral  ratio 
for  a one-layered  35  km  thick  reference  model,  while  the  other 
traces  are  for  different  types  of  two-layered  models.  The 
depth  (in  brackets)  of  the  intermediate  interface  and  the 
P-velocity  above  and  below  it  are  given  on  the  figure. 


Fig.  VII. 3. 3 Illustrating  the  fact  that  a deep  crust  with  high  velocity 
and  a shallower  crust  with  lower  velocity  will  have  the 
peaks  in  the  theoretical  spectral  ratio  at  the  same  fre- 
quencies. In  the  upper  case  the  Moho  depth  is  35  km  and  the 
crustal  P-velocity  is  5.7  km/sec,  while  in  the  lower  case 
the  corresponding  numbers  are  40  km  and  6.6  km/sec. 
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Fig.  VII. 3. 4 
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Illustration  of  the  effect  of  different  window  length  in  the 
calculation  of  the  spectral  ratio.  The  model  used  is  the 
crustal  model  (left  panel)  and  the  whole  crust  and  upper 
mantle  model  of  King  and  Calcagnile  (1976).  The  top  trace  in 
both  cases  is  the  theoretical  ratio,  while  the  other  traces 
give  the  ratio  that  will  be  observed  for  different  lengths 
of  the  Fourier  window.  Window  length  in  seconds  is  given  on 
the  figure.  The  source  function,  the  source-side  crust  response 
the  effect  of  the  mantle  as  well  as  the  effect  of  the  instru- 


ment response  have  been  ignored  in  the  calculations. 
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VII. 4 Variance  of  Teleseismic  P-wave  Amplitudes 

The  large  scatter  between  individual  station  determinations 
for  a given  seismic  event  is  a problem  of  considerable  im- 
portance when  estimating  magnitudes  from  network  data.  Success- 
ful application  of  maximum  likelihood  techniques  (Ringdal, 

1976)  generally  requires  that  the  standard  deviation  c of 
the  world-wide  magnitude  distribution  is  known,  or  at  least 
may  be  assumed  to  lie  within  fairly  narrow  limits.  In  recent 
years,  several  studies  have  been  performed  to  estimate  o 
based  on  world-wide  m^  observations  of  large  events,  most  of 
them  using  WWSSN  data.  The  most  comprehensive  of  these  papers 
(Veith  and  Clawson,  1972)  found  o=0.35  using  a smoothed 
distance-amplitude  relationship.  Evernden  and  Kohler  (1976) , 
also  using  large  events,  have  studied  the  effects  on  a when 
compensating  for  regional  crustal  structure  and  station  site 
geology.  Their  conclusion  is  that,  with  this  procedure,  a value 
of  o=0.21  can  adequately  represent  standard  deviation  on 
a world-wide  basis  for  single  stations,  while  a=0.15  would  be 
adequate  for  observations  by  small  arrays  as  specified 
in  Evernden  (1971) . They  also  argue  on  the  basis  of  these  low 
a values  that  network  magnitude  bias  is  unimportant  in  the  con- 
text of  the  usage  of  an  Mg-m^  discriminant  for  networks  as 
described  in  Evernden  (1971)  . 

An  important  question  is  to  what  extent  world-wide  scattering 
observed  for  large  earthquakes  (typically  m^  = 6.0  and  above) 
is  representative  for  smaller  events,  which  of  course  are  the 
ones  of  most  interest  when  considering  network  bias  problems. 

In  order  to  investigate  this  problem,  it  is  clearly  necessary 
to  employ  another  procedure  than  looking  at  world-wide  magnitudes 
of  a given  event.  This  is  because  the  problem  of  non-detections 
at  individual  stations  becomes  quite  severe  with  existing  earth- 
quake catalogues  for  low-magnitude  events.  Instead,  we  choose 

in  this  paper  the  indirect  approach  of  comparing  m^  values  of 
the  large  arrays  LASA  and  NORSAR  for  events  from  selected 
regions,  assuming  independence  of  their  respective  m^  estimates. 
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The  results  of  this  analysis  are  presented  in  Table  VII. 4.1, 
and  show  that  the  scatter  is  significant  and  appears  to  be 
quite  region-dependent.  For  the  six  geographic  regions  analyzed, 
we  have  estimated  an  average  o=0.30  units  for  the  standard 
deviation  around  the  'true'  magnitude  of  each  array  m^.  At 
individual  regions,  o ranges  from  0.21  (Central  America)  to 
0.40  (Kurile  Islands)  - see  also  Fig.  VII. 4.1.  It  should  be 
noted  here  that  any  systematic  regional  station  bias,  including 
average  signal  attenuation  effects  due  to  local  site  geology, 
have  been  eliminated  prior  to  computing  the  variances.  It 
appears  that  the  scatter  could  be  reduced  if  smaller  regions 
are  considered,’  however,  it  might  be  noted  that  Ringdal  (1974) 
found  a value  of  a as  high  as  0.28  for  an  aftershock  sequence 
from  the  Kurile  Islands  17-30  June  1973,  comparing  NORSAR  and 
LASA  reported  m^  values.  Hence,  a significant  residual  scatter 
remains,  even  when  very  small  source  areas  are  considered. 

Our  a estimates  are  significantly  higher  than  those  of  Evernden 
and  Kohler  (1976)  . This  may  reasonably  be  attributed  to  the 
difference  in  spectral  characteristics  of  large  and  small 
earthquakes,  i.e.,  smaller  earthquakes  tend  to  have  a relatively 
higher  corner  frequency  (Aki,  1967).  In  fact.  Fig.  VII. 4. 2 
shows  the  standard  deviation  of  log  amplitudes  across  the 
NORSAR  array  in  five  frequency  bands  (each  of  bandwidth  0.4  Hz). 
The  data  represent  average  values  of  10  events,  listed  in 
Table  VII.  4. 2.  A clear  increase  in  ampli*  •’ '.t;  scattering  is 
observed  with  increasing  frequency,  with  a standard  deviation 
at  2 Hz  almost  double  that  at  0.6  Hz. 

Hence,  these  results  indicate  that  world-wide  m.  scatter 

b 

increases  with  decreasing  magnitudes,  and  thus  that  magnitude 
bias  studies  based  on  large  events  will  not  necessarily  be 
adequate  at  low  magnitudes. 


F.  Ringdal 
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Fig.  VII. 4.1  Comparison  of  NORSAR  and  LASA  reported  for  50  events 
randomly  selected  from  the  Kurile  Islands  area  (Flinn- 
Engdahl  regions  220-222).  The  straight  line  has  a slope 
of  A=1.0,  and  its  intercept  B (denoting  the  average 
difference  ( NORSAR ) -m^ ( LASA ) ) and  the  orthogonal  standard 

deviation  relative  to  this  line  are  also  specified. 
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N’DARD  DEVIATION  (dB) 


Events  used  for  studying  NORSAR  amplitude  variations. 


TABLE  VII 
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I . INTRODUCTION 

This  note  describes  a new  format  for  the  NORSAR  High  Rate  (HR)  tape, 
capable  of  accepting  data  from  maximum  7 subarrays.  This  format  is 
valid  from  1 October  1976. 


- rj2  - 

Bytes  8-9  Automatic  Subarray  Status 

2 bits  per  subarray,  according  to  the  following 
code : 


Field 


3 


01:  all  sequential  channels  invalid 

10:  all  submultiplex  channels  invalid 

11:  all  channels  invalid  or  subarray  not  defined 

00:  non  of  the  above  (i.e.,  OK). 


2 bits  per  subarray  gives  14  bits  for  7 subarrays, 
bits  14-15  will  be  set  to  11. 

Bytes  10-30  Automatic  Seismometer  Status  4 

24  bits  per  subarray,  according  to  the  following 
code: 

Bit  0-5:  1 bit  for  each  sequential  channel. 

Bit  set  denotes  invalid  channel 
8-23:  1 bit  for  each  submultiplex  channel. 

Bit  set  denotes  invalid  channel. 


7 subarrays  occupy  21  bytes. 

Bytes  31-32  Manual  (on-line)  Subarray  Status 
2 bits  per  subarray,  as  for  Field  3 
(2  unused  bits,  set  to  11) 

Bytes  33-53  Manual  (On-line)  Seismometer  Status 

24  bits  per  subarray,  as  for  Field  4. 

Bytes  54-60  Multisample  Subarray  Indicators 

8 bits  (=1  byte)  per  subarray,  according  to 
the  following  code: 

Bit  on: 

0 line  decommissioned 

1 subarray  in  total  ICW  scan 

2 " being  synchronized 

3 " phone  line  being  tested 

'4-7  spare 

7 subarrays  use  7 bytes. 

Byte  61  Padding  byte,  contains  x'AA' 

Byte  62-323  1st  Logical  Record 


Byte 

324-585 

2nd 

Logical 

Record 

Byte 

586-847 

3rd 

Logical 

Record 

5 


6 


7 


* 
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Byte  848-1109  4th  Logical  Record 
Byte  1110-1371  5th  Logical  Record 


Each  logical  record  constitutes  262  bytes,  and  consists  of  a Status 
Field  of  10  bytes  and  7 pairs  of  Subarray  Data  Blocks  (SDBs) , one 
pair  for  each  subarray,  one  SDB  for  each  sample  period  (each  subarray 
is  sampled  twice  per  decisecond) . The  formats  of  the  Status  Field  and 
of  an  SDB  are  given  below.  See  also  Figure  1. 


The  format  of  the  Status  Field  is  as  follows: 

Bytes  Contents 

0 Length  of  this  status  field  in  bytes,  (currently,  101 

1 Version  number  of  this  status  field.  (Currently,  1) 

2-3  High  rate  sensor  subjected  to  dc  offset  thresholding  as 

follows : 


4-5 

6-8 


9 


Bit 

Description 

0-1 

Reserved  (zero) 

2 

Set  when  sensor  failed  threshold 

3-15 

Sensor  number 

Low  rate 

sensor,  otherwise  the  same 

as  bytes  2 

Discrete 

output  fields  transmitted 

as  follows: 

Bit 

Description 

0-7 

Subarray  number  to  which 

transmitted 

8-15 

DOF  1 

16-23 

DOF  2 

Spare  (zero) 


Any  fields  lacking  entries  shall  be  coded  zeros. 


Each  Subarray  Data  Block  (SDB)  contains  data  for  one  sample  period  from 
orio  subarray.  Two  consecutive  SDBs  are  from  one  and  the  same  subarray. 
The  SDB  format  is  as  follows: 


Fir;  If] 
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By  t*'S^ 
1 0-3 


Byte 

0 


1 


2 

3 


Contents 

Sample  period  indicators,  according  to  the 
following  layout: 

Bi t Meaning  (when  set) 

0 ICW  sync  error  occurred 

1 ICW  poly  error  occurred 

2 ODW  sync  error  occurred 

3 ODW  poly  error  occurred* 

4-7  Spare  (zero) 

0 ODW  received  late  in  previous  sample  period 

and  is  being  treated  as  early  for  current  sample 

1 This  data  sample  is  repeated  from  previous 
sample  period 

2 Data  present  this  sample  period 

3 Redundant  ODW  received 

4-7  Spare  (zero) 

0-3  Spare  (zero) 

4-7  Contains  received  random  data  address  (RDA) 

0-5  Contains  function  route  reflected  from 

subarray  (FRS) 

6-7  Contains  function  select  reflected  from 

subarray  (FS) 


* No  data  present  in  fields  2 and  3.  They  contain  a direct  transcription 

of  the  received  message,  including  polynomial  residual  with  the  received 
sync  pattern  in  field  2,  and  the  rest  of  the  message  in  field  3, 
with  unused  bytes  set  to  zeros.  Bits  0,  1 of  Byte  0,  all  of  Bytes  2-3 
are  zero. 

Field  Bytes  Contents 

2 4-5  Submultiplex  channel  data  value,  defloated  with 

extended  sign,  right- justified  in  bits  0-13. 

Bits  14  and  15  are  zero. 

3 6-17  6 sequential  channel  values,  in  the  same  format 

as  for  Field  2. 

As  indicated  under  the  description  of  Field  1,  the  SDB  format  is  different 
when  an  ODW  polynomial  code  error  occurs. 
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'NR'  RECORD  ID 
BINARY  TIME- 


CONTROL  FIELDS 


HR- BLOCK 


CONTROL  FIELDS 
(62) 

LOGICAL 
RECORD  1 

(262) 

LOGICAL 
RECORD  2 

(262) 

LOGICAL 
RECORD  3 

(262) 

LOGICAL 
RECORD  4 

(262) 

LOGICAL 
RECORD  5 

(262) 

0.5  SECONDS  OF  DATA, 
1372  BYTES 


GENERAL  INDICA- 

_ oZ  bytes 

TORS  1 

AUTOMATIC  SA  STATUS 
_ " _ . SE  " 

MANUAL  SA  STATUS 
” _SE 

MULTISAMPLE 
SA  INDICATORS 


STATUS  FIELD  (10) 

SDBU(18) 

SDB12 

“SDBll  " 

sdb22 

~ sdb31 

SDB32~ 

SDB41 

SDB42 

sdb5i 

SDB52 

SDB61 

SDB62 

sdb7i 

1 

sdb;2 

LOGICAL  RECORD 
0 . 1 SECONDS  OF  DATA , 
262  BYTES 


SAMPLE 

PERIOD 

HTTTBICATOI 


; sub 

3 MUX 

rnsRS" — 


» 8 10  12  14  16  N 

SEQ  1 SEQ  *SEQ  1 SEQ  1 SEQ  'sEQ 
1 ' 2 1 3 ' 4 ! 5 6 


SUBARRAY  DATA  BLOCK 
0.05  SECONDS  OF  DATA 
FROM  ONE  SUBARRAY, 

18  BYTES 


Figure  1:  High  Rate  Block  Layout.  If  4 HR-blocks  comprise  1 physical 

record  on  tape,  then  one  tape  may  contain  3.97  hours  of  data. 
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3.  MODIFICATIONS  IN  TAPE  LABELS  AND  ACR  RECORD 

The  ISRSPS  File  Header  Label  which  starts  the  NORSAR  High  Rate  tape  has 
been  modified  to  indicate  the  format  change  (see  REF  110N,  section  2-1-2). 
Field  7 (characters  39  through  40)  contains  the  EBCDIC  numeric  version 
number  of  the  tape,  which  is  presently  11.  This  number  has  been  changed 
to  21^  to  reflect  the  dramatic  changes  in  the  format  of  the  tape. 


The  parameters  of  the  Array  Configuration  Record  (ACR)  have  been  changed 
as  follows  (see  REF  110N,  section  2-2-1-2) : 


Field 

Bytes 

Previously 

Has  been  changed  to 

1 

2 

22 

7 (number  of  sources) 

4 

6-7 

"02" 

"03"  (source  configuration  level) 

5 

8-9 

4986 

5488  (no.  of  bytes  per  physical  record) 

6 

10-11 

890 

262  (no.  of  bytes  per  logical  record) 

7 

12 

5 

20  (blocking  factor) 

The  source 

data 

description  fields 

need  not  be  changed,  except  for  the 

order  and  the  number  of  such  fields,  within  the  ACR. 


D.  Riaber-Mohn 
F.  Ringdal 
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APPENDIX  B 


FORMAT  OF  NORSAR  LOW  RATE  DATA  TAPE  EFFECTIVE  1 OCTOBER  1976 


niifT 


TAPE 


NORSAR  LOW 
GENERAL  FORMA 
10/15/76 


,OW  RATE  TAPE 


NORSAR 


TAPE  FORMAT 


(Modified  October  1976) 


LOAD  POINT 


VOL-SERIAL  HEADER  RECORD  (80  BYTES) 


HEADER  LABEL  (80  BYTES) 


END  OF  FILE  (TAPE  MARK) 


DATA  RECORD  (1802  BYTES) 


DATA  RECORD  ( 1802  BYTES) 


TRANSATLANTIC  DATA  RECORD 


DATA  RECORD  (1802  BYTES) 


END  OF  FILE  (TAPE  MARK) 


TRAILER  LABEL  (80  BYTES) 

END  OF  FILE  (TAPE  MARK) 
END  OF  FILE  (TAPE  MARK) 


#fR*CfqiNB  BLANK- NOT  riJLKSD 


- 60  - 


NORSAR  LOW  RATE  TAPE 
GENERAL  FORMAT  (2/2) 
10/15/76 


DATA  RECORD*  FORMAT 


•NR' 

i 

BYTES 

0-1 

DATA 

FRAME 

no  1 

BYTES 

2-61 

DATA 

FRAME 

no  2 

BYTES 

62-121 

i 

j 

BYTES 

• • t • 

DATA 

FRAME 

no  28 

BYTES 

1622-1681 

1 

DATA 

FRAME 

no  2 9 

! 

BYTES 

1682-1741 

D 

% 

> 

FRAME 

no  30 

BYTES 

1742-1801 

TOTAL  RECORD  LENGTH 
FRAME  LENGTH 


1802  BYTES 
6 0 BYTES 


NORSAR  LOW  RATE  TAPE 
- 61  - LABEL  FORMATS 

10/1-5/76 


VOL-SERIAL  HEADER 


BYTES 

0-3 

HEADER  IDENTIFIER 

' VCfLl ' 

BYTES 

4-9 

VOL-SERIAL 

e .g . ' Q0Q013 ' 

BYTE 

10 

SECURITY  INDICATOR 

'1' 

BYTES 

11-79 

RESERVED  FOR  EXPANSION 

•b. . .b' 

TOTAL 

LENGTH  IS 

80 

BYTES 

NOTE : 

All  data 

in 

this  record  shoula  be 

interpreted  as 

alphanumeric.  'VOL1'  and  ' 1*  are  fixed.  The  ' 1 ' 
actually  implies  a permanent  retention  cycle.  The 
VOL-SERIAL  represents  the  tape  number. 


HEADER 

LABEL 

BYTES 

0-3 

IDENTIFIER 

' IID’Rl ' 

BYTES 

4-14 

TAPE  TYPE 

’ LOWbRATEbbb  * 

BYTES 

15-20 

LOCATION 

' NDPCbb ' 

BYTES 

21-26 

VOL-SERIAL  (THIS  TAPE) 

e.g.  '000056' 

BYTES 

27-28 

RESERVED 

'bb' 

BYTES 

29-34 

VOL-SERIAL  (PREVIOUS  TAPE) 

e.g. '000055' 

BYTES 

35-36 

LOCATION  CODE 

' 03 ' ' (NDPC ) 

BYTES 

37-38 

TAPE  TYPE  CODE 

'14'  ( LOW  RATE ) 

BYTES 

39-40 

TAPE  FORMAT  VERSION 

' 02' 

BYTES 

41-46 

DATE  = 'bYRDAY' 

e.g.  '1?7Q181 ' 

BYTES 

47-52 

RECYCLE  DATE 

' b99365 

BYTE 

53 

RETENTION  CODE 

'1'  (PERMANENT) 

BYTES 

54-59 

DATE  = 'bYRDAY' 

SAME  AS  BYTES  41-46 

BYTES 

60-65 

PROCESSING  SYSTEM 

' ISRSPS ' 

BYTES 

66-79 

COMMENTS 

' NORSARbLPbDATA ' 

TOTAL  LENGTH  IS  80  BYTES 

NOTE:  All  data  in  this  record  should  be  interpreted  as 

alphanumeric  characters. 


TRAILER 

LABEL 

BYTES 

0-3 

IDENTIFIER 

' EOF1 ' 

BYTES 

4-28 

RESERVED 

'bb. . .bb' 

BYTES 

29-34 

VOL-SERIAL  (NEXT  TAPE) 

e.g.  '000057' 

BYTES 

35-79 

RESERVED 

' bb. . . .bb 

- 62  - 


wOKSAR  LOW  RATE  TAPE 
DATA  FRAME  FORMAT 
10/15/76 


DATA  FRAME  FORMAT 
LENGTH 


11 

POSITION 

5 BYTES 

RAW  TIMING  WORD 

0-4 

1 BYTE 

RAW  SITE  ERROR  TAB!..: 

5 

4 EYTES 

RAW  SEIS  FUNCTION  TABLE 

6-9 

42  3YTES 

RAW  LONG  PERIOD  DAS A 

10-51 

8 BYTES 

FRAME  FLAG  BYTES 

52-59 

?OTAL  LENGTH  IS  60  BYTES 


RAW  TIMING  WORD 

FRAME 

BYTE 


(5  BYTES 


NOTE : 


i i ‘ ' 


V) 

-p 

u 

K 

o; 

b ^ 

0 

-P 

P 

(/) 

o 

4-1  W 
3 -P 

C 

O 

a: 

>i*H 

3 

C 

— ' ,r"> 

£ 

c 

C -H 

o 

t: 

0 

0) 

O £ 

•H 

o 

■H  C 

<u 

o 

Tj  P 

4-> 

.c  s 

U 

-P 

S 3 

05 

-p 

Digits  in  timing  word  are  single  HEX  digit-BCD 
e.g.  00102  = 21q  = 2 16 


10012  " 910 


= 

16 


HEX  characters  A thru  F are  illegal 


Zeros 


RAW  SITE  ERROR  TABLE  ( 1 BYTE,.  5) 
FRAME  . 


FRAME  V 
3YTE  J 


0 

T 

1 

— r 

2 

1 1 * 

3 4 5 6 

"1 

7 

1 

CM 

n 

o 

o 

0 

0 

0 

p 

c 

c 

c 

c 

P 

a) 

a) 

N 

mJ 

U) 

• • • w 

cn 

£-• 

f4 

Eh 

3 

M 

M 

►H 

M 

C 

Oj 

05 

o. 

05 

t> 

NOTE:  One  bit . represents  one  site. 

Bit  Value  = 0 : normal  Bit  Value  = 1 ; abnormal 

This  bit  configuration  is  to  be  referred  to  as  the  SI" 
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NORSAR  LOW  RATE  TAPE 
DATA  FRAME  FORMAT 
10/15/76 


RAW  SEISMOMETER  FUNCTION  TABLE  BYTES  6-9) 


FRAME 

BYTE 


NOTE:  Four  bits  represent  one  site-2  site;;  are  represented  in  one  byte. 

Bit  0 or  4 corresponds  to  the  Vertir  1 component  (V) 

Bit  1 or  5 corresponds  to  the  North-,  outii  component  (NS) 

Bit  2 or  6 corresponds  to  the  East-N.  st  component  (EW) 

Bit  3 or  7 unused 

Bit  Value  = 0 : normal  site/ccmpr  ent 

Bit  Value  = 1 : abnormal  site/compt „ent 


9 _ 1.0  ^^11 


/ V 

■ia  aW 

1 

.4  * 1UI 

14, 

I H J ■■ r 

rH 

<N 

ro 

r-' 

o 

o 

O 

o 

c 

c 

C 

c 

TJ 

a) 

H 

w 

cu 

w 

0J 

9-4 

(r< 

£-< 

E-i 

3 

H 

H 

M . 

M 

C 

w 

Vi 

w 

cn 

D 

NOTE:  This  configuration  of  bits  will  be  i • noted  SFT  format. 


RAW  LONG  PERIOD  DATA  (42  BYTES  10-51) 


NOTE:  Long  Period  Data  is  recorded  on  tape  as  integer  halfwords. 

There  are  tv:o  bytes  per  sample  point.  The  LP  Data  contains 
7 sites  with  three  components  per  site  (the  components 
vary  most  rapidly) . 


The  Long  Period  Data  values  differs  from  those  recorded  on 
the  NDPC  data  tapes  by  not  having  the  low  order  two  bits 
used  for  error  checking.  Thus  before  recording  on  the  Low 
Rate  tape  takes  place  the  data  values  received  from  SPS  are 
d ivided  by  4 . 
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NORSAR  LOW  RATE  TAPE 

DATA  FRAME  FORMAT 
STATUS  FRAME  FORMAT 
10/15/76 


FRAME  FLAG  BYTES 


(8  BYTES  52-59) 


These  eight  bytes  are  unused  and  will  always  contain  zeros. 


i 

: 


r 
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NORSAR  beam  set  valid  from  1 October  1976.  The  first  part  of 
the  table  is  the  coherent  beam  set  (beam  set  451,  180  beams), 
while  the  second  part  is  the  incoherent  beams  (beam  set  271, 

64  beams) . UX,UY  given  in  the  table  corresponds  to  the 

slowness-space  location  where  the  phase  given  (PHASE)  is 

expected  to  be  observed  for  an  event  taking  place  at  the 

latitude  (LAT)  and  the  longitude  (LON)  listed.  Region  number 

and  name  correspond  to  the  Flinn  and  Engdahl  (1965)  regionalization. 

H.  Bungum 
F.  Ringdal 
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I 


NOMSAR  ARRAY  BEAM  $tT  <.51  UCTUa^R  l.l  9/c 


BEAM  NC 

UX  ( S/KM  1 

UY  IS/KMI 

PmA  jt 

l A 1 

L UN 

Mt  UN 

NljM|*f  8 ANO  NAME 

1 

C. 008990* 

-0.091 390. 

p 

o 3N 

7H 

6**  . 

NtJK  1 M <lf  SVALUAKD 

2 

-0.0638110 

-C. 08737 CO 

P 

7 3N 

556 

6*  H 

NUVAYA  Z t Ml  YA 

J 

c.oaa  360 : 

-G.06893CC 

P 

72  4 

7A  m 

t*|2 

«A|  F IN  ISLAND  REGION 

A 

C.C1 72A0C 

-G . 0o5  90CC 

P 

6 7n 

1 376 

67/ 

N.  Yt/KUN  TERR.  • CA9ACA 

5 

0 • Go  9 7 601 

-C.C65U00C. 

P 

6 3N 

151m 

1 

l.l  NT  HAL  ALASKA 

6 

-.01177  'C 

-0.C6  IHlOi 

P 

6 1 N 

1A8W 

2 

SOUTMt-  RN  ALASKA 

7 

C.OC^SttOC 

-0. J6*5i99 

P 

65N 

1 65  W 

bio 

ALASKA 

o 

O.C 18A8CC 

-C  . 0 59980’J 

P 

59N 

137m 

19 

SOUTHEASTERN  ALASKA 

9 

0.01 3300C 

-0.061  750«; 

P 

6^N 

156M 

2 

SlIU  T m RN  ALASKA 

10 

O.OC5A9Jr 

-0.0602300 

P 

5 7N 

158H 

12 

ALASKA  PENINSULA 

11 

J.0o93C0i 

— J.  0596800 

P 

5 7N 

1 52  M 

1 3 

KOu|AK  ISLANO  86  6 f ' *N 

12 

C. 0019530 

-g.(  5757CC 

P 

55N 

l 6 3 ta 

1 

UNIMAK  island  region 

li 

-0.G3173OC 

-0.0581200 

P 

6 2N 

15  36 

671 

6 AS  T C 8 i SIBERIA 

1A 

-0.00  30  30r 

-0.0565900 

P 

52N 

171m 

9 

FUX  ISLANDS.  ALEUTIANS 

15 
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